Transport in disordered graphene nanoribbons 
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We study electronic transport in graphene nanoribbons with rough edges. We first consider a 
model of weak disorder that corresponds to an armchair ribbon whose width randomly changes by 
a single unit cell size. We find that in this case, the low-temperature conductivity is governed by 
an effective one- dimensional hopping between segments of distinct band structure. We then provide 
numerical evidence and qualitative arguments that similar behavior also occurs in the limit of strong 
uncorrelated boundary disorder. 
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I. INTRODUCTION 

Since the invention of the method for production 
graphene 1 , many creative ideas for physical effects and 
devices have been put forthS. Whereas early papers em- 
phasized unusual electron properties of graphene as com- 
pared with ordinary metallic and semiconductor materi- 
als, it had been soon realized that graphene is a promising 
material for implementation of previously known physical 
devices with considerably improved characteristics. One 
of possible applications would be in semiconductor tech- 
nology: excellent mechanical properties, easily tunable 
electron concentration, zero nuclear spin, and simple pro- 
duction are among the advantages that make graphene 
a much sought after material. However, the drawback 
preventing the use of graphene for semiconducting ap- 
plications is exactly what is usually considered to be its 
main feature — the absense of a gap in the spectrum. In 
the absence of the gap, it is impossible to make even the 
simplest conventional electronic devices. For instance, 
a graphene p-n junction does not rectify current, even 
though it has some other interesting properties due to 
the Klein tunneling^. Similarly, hybrid graphene - nor- 
mal metal systems are conducting for arbitrary gate volt- 
age applied to graphene^. The only known way to open a 
gap in monolayer graphene is to use confined geometries 
- graphene quantum dots^ and graphene nanoribbons 
(GNR). 

The electronic structure of ideal GNR is theoretically 
well established. It is very sensitive to the ribbon geome- 
try, i.e. orientation relative to the crystal axes and their 
exact widtbJii£&2. Within a tight-binding model with 
only nearest neighbor hopping, GNR with zig-zag edges 
have flat near-zero-energy bands of extended edge states, 
while ribbons with armchair edges, depending on the pre- 
cise width, can be either metallic or semiconducting with 
the gap inversely proportional to the GNR width. Nu- 
merical studies^ show that passivation of the edges of 
ideal GNR — chemical bonding of edge carbon atoms 
with hydrogen — may open a small energy gap that does 
not scale with the GNR width. 

Recently, first experimental observations of transport 
in GNR have been reportediiJ^I^. GNR with 



widths in the range from about 10 nm to 100 nm and 
lengths in the micrometer range have been studied. The 
fabrication procedure does not yet allow to control the 
GNR width with atomic precision (although chemical 
fabrication-^ may eventually yield controlled edge fab- 
rication). As a result the edges are disordered on the 
atomic length scale, as well as show longer-range width 
variation of a few nanometers. For narrow enough rib- 
bons 50 nm) an unambiguous signature of the geomet- 
ric gap Eg scaling with the inverse average ribbon width 
has been extracted from the gate voltage and tempera- 
ture dependencies of conductivit y 11 ' 12 . In particular, in 
Ref.flll in a broad range of temperatures, T, conductivity 
scales as e~ Eg ^ T . The measured gap is a smooth func- 
tion of the ribbon width, and is insensitive to the GNR 
orientation relative to the crystal axes. Also, 1// current 
noise has been observed at low frequencies, / < 100 Hz, 
with the intensity proportional to the GNR widthii. 

These experimental results are definitely inconsistent 
with the theory for ideal GNR that predicts different be- 
havior depending on the orientation, typically with many 
low-energy states. The observed effects are clearly due 
to disorder. Indeed, it is natural to expect that any dis- 
order, bulk or boundary, should lead to Anderson local- 
ization and open a transport gap; however, one would 
expect that this gap should be defined by the strength of 
disorder, rather than by the GNR width. 

In this work we provide a qualitative resolution to this 
apparent puzzle by showing that electronic properties of 
the disordered GNR are indeed very different from the 
clean GNR. We demonstrate that for the states near the 
middle of the band, edge disorder leads to segmentation 
of the wavefunctions into blocks of length of the order 
of GNR width. Thus, at low temperatures, the sys- 
tem maps onto an effective one-dimensional (ID) hop- 
ping insulator—. We illustrate this behavior first with a 
model where disorder is introduced through slowly fluc- 
tuating ribbon width, which allows more direct numerical 
and analytical analysis (Section [TTJ) . Then we generalize 
the results to the experimentally relevant case of strong 
disorder (Section IIIip . Discussion and conclusions are 
presented in Section HVl). 
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II. WEAK DISORDER 




FIG. 1: Structure of the electronic wavefunctions in a "weakly 
disordered" armchair GNR. Periodic boundary conditions are 
applied in the horizontal direction. In our convention, the 
narrow segments have width W = 4\/3, and in the limit of 
infinite length would have no band gap; the surrounding re- 
gions (width W = 5\/3) for an infinitely long ribbon would 
have a gap E g = 2 x 0.169£ around zero energy (the middle 
of the GNR 7r band). The radii of the circles are propor- 
tional to the site amplitudes of the wavefunction, with the 
color representing sign. The top plot corresponds to a the 
lowest energy state inside the gap, E = — 0.096t, spatially lo- 
calized in the left (longer) metallic segment; the middle - to 
the low energy state in the right (shorter) metallic segment, 
E — — 0.13H; the bottom is a delocalized state well outside 
the gap, E = —0A97t. Note the abrupt change in the local- 
ized wavefunctions' amplitude at the "interface," and rather 
uniform amplitude across the ribbon. 

Let us consider an armchair GNR. An ideal ribbon 
of the width W, measured in units of minimal carbon- 
carbon distance, a g , is metallic (no gap) for W = (3N + 
1)V3, and semiconducting (with gap Eg ~ t/W) for 
W = 3VSN and W = (3N + 2)%/3, Ref. Here, t is the 
graphene nearest-neighbor hopping matrix element (we 
neglect the next-nearest neighbor hopping which causes 
slight particle-hole asymmetry), and N is an integer. 
Weak disorder can be introduced as geometric fluctua- 
tions of the ribbon width, such that the "disordered" 
ribbon is comprised of ideal segments of random length of 
order L, with width changing from segment to segment. 
We assume that L > W . An example of a "disordered" 
configuration of this kind is shown in Fig. [TJ While this 
situation has not been yet realized experimentally, it has 
the advantage that its analysis is straightforward, and, 
as we will argue, the behavior is related to the experi- 
mentally relevant case of strong disorder. If the length 
of each segment is longer than its width, to the low- 
est order, one can consider individual band structure of 
each segment separately. Depending on the width, some 
of the segments are nearly metallic, with the finite size 
gap of about t/L, while others are "insulating" with the 
gap t/W. In Fig. [1] we show the results of numerical 
diagonalization of the tight-binding graphene Hamilto- 
nian — representative wavefunctions in different regions. 



Note that although the segmentation is caused by the 
surface defect (change of the width by just one row of 
atoms!), the wavefunctions show high degree of unifor- 
mity across the ribbon, and rather sharp confinement to 
the respective regions along the ribbon. Thus, at low en- 
ergies (\E\ < t/W) it is natural to represent the system 
by a one-dimensional hopping model, 

H = J2^fcT+J2^f^ + h-c. (1) 

ia ia,j@ 

Here, the operator c"^ (cf ) creates (destroys) electron in 
the metallic "grain" i in the orbital a. 

To complete the formulation of the effective model 
Eq. ([T]) we need to determine the distributions of the on- 
site energies ef and inter-site hopping matrix elements 
tf.® . For simplicity we assume that the average length of 
the segments, both insulating and metallic, is the same, 
L av . The low-energy spectrum in the metallic segments 
follows from the Dirac dispersion of the infinite metallic 
armchair GNR2, e = c\k\, where c = 3ta g /2 and k is the 
momentum along GNR. The levels in a given metallic 
segment of length L are therefore approximately equidis- 
tant, with the average level spacing ~ t/L. In Figure[2]we 
show the result of a tight-binding calculation for the low- 
est energy state as a function of the length of a metallic 
segment embedded in the insulating GNR. Indeed we find 
that the energy scales approximately as 1/ L. Even better 
fit is obtained by using the form 1/(L + L\y) which takes 
into account the leakage of the wavefunction from the 
metallic regions into surrounding insulating ones. Note 
that L\v ~ W. If the lengths Li for all grains were equal, 
the level structures in all grains would be identical (apart 
from the small splitting caused by inter-grain tunneling) . 
However, for a distribution of lengths, the energy levels 
in different grains are likely to be out of registry by the 
amount ~ t/ L av . 

We now evaluate the tunneling matrix elements ttj be- 
tween low-energy states in metallic segments. Tunneling 
occurs through the intermediate states in the insulat- 
ing regions. The states just outside the gap are par- 
ticularly important for tunneling. Near the gap edge 
the dispersion is quadratic, e = y/c 2 k 2 + (E g /2) 2 w 
Eg/2 + c 2 k 2 / Eg. This corresponds to the effective mass 
in the insulating regions in* ~ (Wa^t) -1 . The tunnel- 
ing amplitude through a barrier of height E g and length 
D can be estimated using the WKB approximation as 
e ~ aD / w j where a is a numerical coefficient of order 1. 
We have also verified this by a direct tight-binding cal- 
culation of the tunnel splitting of energy levels in two 
identical metallic segments separated by an insulating 
segment of variable length. 

From the distributions of ef and tff, it follows that 
for L av > W, the level spacing in the metallic grains is 
larger than the tunneling amplitude between the neigh- 
bors, making it impossible to have metallic, i.e. band, 
conduction. The system is a one-dimensional example of 
a simple impurity band insulator, a standard model used 
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FIG. 2: The lowest energy state in a GNR as a function of 
length L of a "metallic" constriction surrounded by an "in- 
sulator." The configuration is similar to the one in Fig. [T] 
The total length of the ribbon used in simulation is 120a 9 
(with periodic boundary conditions along horizontal axis). 
The results presented are for two GNR widths, W = 5^3 
and W = llV3, with W = 4^3 and W = 10-^/3, respectively, 
in the metallic regions. In the inset we fit the energy to the 



form E oc (L + L\ 



The offset Lw appears due to the 



leakage of the wavefunction from the metallic regions into the 
insulating ones. As expected (see text), Lw ~ W. 



to describe lightly doped compensated semiconductors 17 . 
The finite-temperature conductivity of such insulator can 
be evaluated by standard techniques 16 , 



-anL a v/W-t/(nL a .vT) 



(2) 



where T is the temperature, and n is the length of the 
optimal jump. By minimizing the exponent, we find that 
n opt = y/tWj (aL^ v T). Hence there is a crossover from 
the nearest neighbor (n opt — 1, NNH) to variable range 
hopping (n^ > 1> VRH) at temperature T* ~ tW/L 2 av , 



,-2y/ a t/(WT) 
D ~aL av /W~t/(L av T) 



for T <T* 

for T* <T < t/L a 



(3) 



Note that T* < t/W, and thus both behaviors are pos- 
sible within our model. At temperatures higher that 
t/L av multiple states in the metallic regions have to be 
included. We do not consider here other regimes of one- 
dimensional hopping^ that can become relevant at very 
low temperatures. 



III. STRONG DISORDER 

We now turn to the strong disorder case, when the 
boundary is randomized at the atomic scale (this models 
the situation when some atoms are cut out or replaced by 
other atoms, e.g. oxygens, in the process of fabrication). 
Yet, we assume that the relative variation of the ribbon 
width introduced by disorder is small. This is different 
from the near-granular case considered in Ref. [T^ . 

An example of a "strongly" disordered configuration 
that corresponds to small relative variation in the rib- 



bon width is shown in Fig. [3] We chose a perfect zig- 
zag nanoribbon as the reference structure. Ideal zig-zag 
nanoribbons are always metallic, for any ribbon width, 
due to the presence of the edge states 6 .. We observe that 
edge disorder (here generated by eliminating at random 
half of the sites along the edges) leads to the wavefunc- 
tion localization. However, since disorder is now short- 
correlated, the wavefunctions no longer have a typical 
extent along the ribbon, but rather can be either more 
or less localized. We find numerically that the wavefunc- 
tions corresponding to the low-energy states (\E\ < t/W) 
that are highly localized along the direction of the bound- 
ary (L -C W, e.g. Fig. [3^) also do not penetrate deep 
inside the ribbon, having large amplitude only near the 
surface. On the other hand, states that are more ex- 
tended along the ribbon also penetrate deeper into the 
bulk. This effect can be traced back to the behavior of 
the edge states in zig-zag GNR - the wavevector along 
the ribbon for these states is approximately equal to their 
exponential decay length into the bulk. In effect, in the 
absence of next-nearest-neighbor hopping, one can have 
states with very low energy, ~ te~ w l a , localized over the 
distance of about single unit cell near the ribbon edge. 
The number of such states within a segment of length W 
can be easily estimated to be \K — K'\/(2n/aW) ~ W, 
that is, each boundary atom can support one highly lo- 
calized low energy state. Their spectrum which can be 
derived from the dispersion of the ideal ribbon surface 
states is2 E ~ texp(-kW) for kW 3> 1, where k is the 
wave vector deviation from the Dirac point. This spec- 
trum provides the density of states which behaves as E~ x 
at E<.t/W. 

This extremely high low-energy density of states is ob- 
viously an artefact of our model and disappears if the 
hopping between next-nearest neighbors (NNN) on the 
graphene lattice is taken into account. It is character- 
ized by the overlap intergral t' which in graphene ap- 
proximately equals to 0.2t. In the following, we assume 
t 1 > t/W which is the case for all graphene nanoribbons 
studied in the experiments. If the NNN overlap is taken 
into account, the edge states are hybridized and form 
a band of the width ~ t'. Thus, the (one-dimensional) 
density of states in the gap \E\ < t/W is approximately 
constant and equal to one state per surface atom per t' . 
Since the localization length is governed by the energy 
distance to the next subband, it is nearly independent 
of t' and still is approximately equal to W. We there- 
fore coarse-grain over the ribbon elements of size W xW 
to obtain the low-energy level spacing within a coarse- 
grained element AE ~ t'/W. Now, we can formulate 
the variable range hopping conductivity between these 
elements as 



-cm-t' /(nWT) 



Optimizing over n, we find two regimes, 



,-2y/at'/WT 
,-t'/{WT) 



for T <T C 

for T C <T < t/W 



(4) 



(5) 
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with the crossover temperature T c — t'/W < 
t/W. Again, higher-temperature (top line) and lower- 
temperature (bottom line) ranges correspond to VRH 
and NNH, respectively. 

On the experimental side, Chen et alM find that, e.g. 
in 20 nm wide GNR, which has a gap 28 mV, at relatively 
high temperatures, between 50 K and 100 K, transport 
is activated, a cx e~ Eg / T . The data is lacking at interme- 
diate temperatures; however, the single low-temperature 
data point at 4 K shows conductivity much higher than 
would be expected from simple activated hopping. This 
may reflect a crossover from the nearest to variable range 
hopping. Our estimate for the crossover temperature is 
about 5 Kelvin, which is consistent with the experimental 
results. More detailed experimental data in the interme- 
diate temperature regime should allow direct test of our 
predictions. It has also been found that the experimen- 
tal size of the gap is smaller than t/W, which may be 
consistent with our value t'/W. 




FIG. 3: Structure of the electronic wavefunctions in a strongly 
disordered zig-zag GNR. The disorder is generated by ran- 
domly eliminating half of carbon atoms at the edges of GNR. 
Periodic boundary conditions are applied in the horizon- 
tal direction. The energies of the states, from top down, 
E = -0.071t, -0.089*, -0.255t. Note how the confinement 
length increases away from the center of the band (E — 0). 
The typical confinement length at the energies inside the gap 
is of the order of the ribbon width. 



IV. DISCUSSION AND CONCLUSIONS 

The alternative explanation of experimental results 
was proposed in a model respresenting a GNR with very 
strong interface disorder as a chain of quantum dots, host- 
ing localized electron states^. It was suggested by the 
authors that the experimentally observed gap has to do 



with charging energy of the "quantum dots." Indeed, re- 
cent experiment al 15 indicate that there are local charg- 
ing centers in graphene nanoribbons. Compared to the 
case studied in Rcf. both cases that we considered 
above correspond to at most mild disorder. We therefore 
address now how the Coulomb interaction will affect our 
results. 

The Coulomb interaction can modify low-temperature 
hopping conductivity. Interaction leads to opening of the 
soft Coulomb gap around the Fermi surface, which enters 
as an energy cost inversely proportional to the length of 
the hopii. Thus in the presence of the Coulomb interac- 
tion, the expression for conductivity has to be modified 
as 

a ~ e -an-t' /(nWT)-e 2 /(ena g WT) _ /g) 

Since e 2 /a g t ~ 1 in graphene, the Coulomb cost will be- 
come relevant if the dielectric constant of the embedding 
medium e is smaller than t/t' ~ 5. While the functional 
form of conductivity in this case remans the same as in 
Eq. ([3]), the energy scale that defines the gap is different, 
t'/W -f e 2 /(ea g W). 

Thus, for freely suspended graphene the transport is 
in fact expected to be dominated by the soft Coulomb 
blockade. On the other hand, placing graphene in the 
vicinity of high-e medium or metallic gate would reduce 
the Coulomb interaction strength and ranged, leading to 
crossover to Mott's VRH. 

Finally, we note that the 1/f noise observed by Chen 
et alJ^ may also be consistent with the scenario pre- 
sented here, that is, it may be intrinsic, rather than 
caused by the charge fluctuations in the substrate, as 
was suggested in Ref. [ill Due to the presence of 
an exponentially broad distribution of the tunneling 
rates in the hopping transport, the experimentally ob- 
served Hooge relational between the low-frequency cur- 
rent noise and the DC current, I 2 /I 2 = A(uj,T)/u>, 
can be naturally expected 22 . Straightforward applica- 
tion of the Shklovskii's argument 22 to one dimension 
leads to Hooge's parameter A cx exp(— BT) in the low- 
temperature (VRH) regime, and approximately constant 
A at high temperatures (NNH). Whether the 1/f noise 
is indeed intrinsic can be tested by varying the substrate 
properties, or performing measurement on a suspended 
GNR 2 !. 
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